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We present a measurement of asymmetries in the production of and A~ 
baryons in 500 GeV/c 7r~-nucleon interactions from the E791 experiment 
at Fermilab. The asymmetries were measured as functions of Feynman x 
(xp) and transverse momentum squared using a sample of 1 819=b 62 
A^'s observed in the decay channel A+ — > pK'ir'^ . We observe more A+ 
than A~ baryons, with an asymmetry of (12.7 =b 3.4 it 1.3) % independent 
of Xp and p| in our kinematical range (-0.1 <xp < 0.6 and 0.0 < < 
8.0 (GeV/c)^). This A^ asymmetry measurement is the first with data in 
both the positive and negative Xp regions. 



Several studies [Q have reported an enhancement in the production rate of 
charm particles having a valence quark or diquark in common with the inci- 
dent particles, relative to charge-conjugate particles which have fewer or no 
common valence quarks. This effect, known as the leading-particle effect, can- 
not be accounted for by next-to-leading-order perturbative QCD (pQCD) cal- 
culations followed by Peterson fragmentation of the produced c and c quarks. 
Thus, the effect is attributed to additional hadronization effects or non-pQCD 
contributions to charm production. 

In this work we report on the — production asymmetry measured using 
a sample containing 1 819 ± 62 fully-reconstructed — > pK~i[^ and charge 
conjugate (c.c.) decays from the Fermilab E791 experiment. The data are from 
500 GeV/c 7r~-nucleon interactions. This measurement represents a significant 
increase in statistics from previously published results, and it also explores for 
the first time the production asymmetry at negative values of Xp. This feature 
of our measurement allows us to look both for diquark effects which should 
be present in the negative Xp region (since the A+ has two valence quarks 
in common with the target particles whereas the Ar has none) and for lead- 
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ing particle effects in the positive Xp region (in which both the particle and 
antiparticle share a single valence quark with the incident pions). The pro- 
duction asymmetries are measured as functions of Xp and p"^ in the kinematic 
ranges -0.1 < < 0.6 and < < 8 (GeV/c)^. 

Varied predictions about what asymmetries should be observed in the E791 
data are made by three models which address leading particle effects: the 
string fragmentation (SF) model 0], the intrinsic charm (IC) model 0, and 
the two-component recombination (2CR) model f^. 

In the string fragmentation model implemented in the Lund P YTHIA / Jetset 
package P, the c and c quarks produced by leading order QCD diagrams are 
connected by color strings to valence quarks or diquarks in the beam or target 
particles. Hadronization proceeds through fragmentation of the ends of the 
color string and recombination with valence quarks from the initial particles. 
The model predicts a negative asymmetry (more A~ than A^) for > 0, a 
positive asymmetry for Xp < 0, and a fiat dependence on p^.. 

For the intrinsic charm model, fluctuation of one or the other of the incident 
particles to a Fock state containing cc quarks {\uudcc) or \uddcc) for nucle- 
ons and \udcc) for vr") favors production of particles with valence quarks in 
common with the Fock state via a coalescence (recombination) mechanism. 
Hence, this model predicts no asymmetry in the > region and a positive 
asymmetry for < 0, increasing as one moves toward more negative Xp. 

In the two-component recombination model, hadronization takes place through 
fragmentation of perturbatively produced c and c quarks, and the recombi- 
nation of valence and c (c) quarks produced in the initial collision. In the 
Xp < region, this model predicts a positive asymmetry growing strongly as 
Xp becomes more negative, due to target-diquark effects. Again, there is no 
asymmetry in the pion fragmentation region (x^ > 0) since both A+ and A~ 
are singly-leading particles. 

The main difference between the IC and 2CR models is that the former pre- 
dicts very little asymmetry in the E791 negative Xp range because the IC 
distribution peaks at about 0.6 of the initial particle momentum. Because of 
diquark effects, the 2CR model can account for larger asymmetries in this 
region. 

Neither the IC nor the 2CR model has been used to make a prediction regard- 
ing the variation of asymmetry with respect to p"^. Whatever asymmetries do 
exist, however, are expected to decrease as p^, increases, since the density of 
beam and target valence quarks available for recombination decreases with 
increasing p^.. 

Experiment E791 recorded data from 500 GeV/c vr^ interactions in five thin 
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foils (one platinum and four diamond) with center-to-center separations of 
about 15 mm. Each foil was approximately 0.4% of a pion interaction length 
thick (0.5 mm for the platinum foil and 1.6 mm for the carbon foils). The E791 
spectrometer in the Tagged Photon Laboratory was a large- acceptance, 
two-magnet spectrometer with eight planes of multiwire proportional cham- 
bers (MWPC) and six planes of silicon microstrip detectors (SMD) placed 
upstream of the target for beam tracking. Downstream of the target were 17 
planes of SMD's for track and vertex reconstruction, 35 drift chamber planes, 
two MWPC's, two multicell threshold Cerenkov counters, electromagnetic and 
hadronic calorimeters, and a muon detector. The experiment recorded 2 x 10^° 
interactions using an open transverse-energy trigger and a fast data acquisition 
system 0. 

A sample of candidate pK^n^ and c.c. decays was selected using similar 

criteria to those reported in Ref. Since the mean decay length for A^'s was 
several mm, in most cases they decayed in air gaps between the target foils 
and before entering the silicon vertex detectors. We selected all combinations 
of three tracks which had Cerenkov counter identification probabilities consis- 
tent with the pKtt hypothesis and which had a pKn invariant mass between 
2.15 and 2.45 GeV/c^. We then required that the distance between the recon- 
structed candidate A^ decay vertex and the production vertex be at least five 
times the rms uncertainty in that distance. To further reduce background, we 
required the Ac to decay a distance greater than da^ from the nearest mate- 
rial, where is the uncertainty on the vertex z position. The vertex position 
is also required to lie between one and four A^ lifetimes from the primary 
(interaction) vertex. The A^ momentum vector, reconstructed from its decay 
products, was required to pass within 50 fim of the production vertex, and the 
reconstructed production and decay vertices were each required to have an ac- 
ceptable x^- We also required that no more than one of the A^-candidate decay 
tracks be consistent with coming from the primary interaction point. The sum 
of the of the secondary tracks with respect to the flight path of the recon- 
structed Ac was required to be greater than 0.35 (GeV/c)^ if —0.1 < Xp < 0.3, 
and greater than 0.32 (GeV/c)^ if 0.3 < a;^ < 0.6. 

Events were kept if the A^ candidate had —0.1 < Xp < 0.6 and < 
8 (GeV/c)^. The pKtt invariant-mass plots for the final A+ and Aj samples 
are shown in Figs. |I](a) and (b) respectively for < 0, and in Figs, ^(c) and 
(d) for Xp > 0. Fits to a Gaussian signal and a quadratic background yield 
1 025 ± 45 A+ pK-'K+ and 794 ± 42 A; ^ pi^T+Tr" in the entire Xp and pi 



range. 



For each bin of Xp or p^, we define an asymmetry parameter A as 



A = 



N -N/r 



e 

r = - 
e 



(1) 



N + N/r 
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where (A^) is the number of (Aj) produced in the bin and the efficiency 
e (e) is the product of the geometrical acceptance and reconstruction efficiency 
for the bin. 

Values for (A^) were obtained from ffis to pK-K invariant-mass distribu- 
tions for events within specific Xp and ranges. In each case, well defined 
Ac signals were evident. Efficiencies (e, e) were calculated using a sample of 
7 X 10'' Monte Carlo (MC) events produced with the Pythia/Jetset event 
generator These events were projected through a detailed simulation of the 
E791 detector and then reconstructed with the algorithms used for the data. 
The final reconstructed MC sample was approximately ten times the size of 
our data sample. In the simulation of the detector, special care was taken to 
accurately simulate the reduced chamber efficiency seen in the experiment in 
a small region of the tracking chambers nearest the beam. The behavior of the 
apparatus and details of the reconstruction code changed during the data tak- 
ing and data processing periods, respectively. To account for these effects, we 
generated the final MC sample in subsets mirroring these behaviors and com- 
bined the samples using their known fractional contributions to the data set. 
Good agreement between MC and data samples in a variety of kinematic vari- 
ables and resolutions was achieved. The efficiency ratios obtained are shown 
in Fig. |. 

Acceptance-corrected asymmetries in the range —0.1 < Xp < 0.6 integrated 
over all < 8 (GeV/c)^, and in the < 8 (GeV/c)^ range integrated from 
-0.1 to 0.6 shown in Fig. |^ and listed in Table |1|. The statistical errors 

given include those due to the number of observed events and the number 
of MC events accepted. Typically, only for the highest Xp and p| bins is the 
error dominated by the MC uncertainty. Also shown in Fig. |^ are predictions 
from the default Pythia/Jetset . The asymmetries predicted by Pythia 
are negative for x^- > and in all the p\ range studied, while our data exhibit 
a positive asymmetry throughout the kinematical region studied. However, 
the data do not exclude a rise in the small negative Xp region. 

Sources of systematic uncertainties were checked. Among them we checked the 
effect of changing the parametrization of the signal and background shape, the 
effect of varying our principal selection criteria and the effect of a 2.5 % K~ 
contamination in the beam. 

The most significant effect comes from the parametrization of the background 
shape and the variation of our principal selection criteria. The effect of the 
contamination in the beam, which might produce a negative asymmetry 
in the Xp > Q region, was found to be negligible. 

We also checked for contamination due to D"^ and Df decaying into i^'vrvr 
and KKti modes. We found that, by restricting our sample to have lifetimes 



5 



between one and four lifetimes, we reduced any possible contamination to 
negligible levels. 



The SELEX collaboration has measured an asymmetry in the Xp > 0.2 region 
using different incident beam particles. Their preliminary results for the vr" 
beam indicate an asymmetry of A = (25 ± 15)% 0. The ACCMOR Collabo- 
ration has also measured the A+ — A^ asymmetry over the < Xp < 0.8 
region in 230 GeV/c 7r~-Cu interactions and found A = (0.5 ± 7.9)%, indicat- 
ing no asymmetry, although with large uncertainty. Our measurement, when 
averaged over the < Xp < 0.6 region, is A = (12.3 ± 3.7 ± 1.6)%. No in- 
consistencies are observed in these results. The asymmetries measured by the 
three experiments and predictions of theoretical models are shown in Table |^. 



In the forward Xp region, leading particle effects are expected to cancel since, 
in this region, both particle and anti-particle share one valence quark with inci- 
dent pions. The observed constant asymmetry could be due to different energy 
thresholds for particle and antiparticle in associated production of charmed 
mesons and baryons. A similar effect was observed in the hyperon production 



asymmetries measured in this experiment [|Tl|]. Our data cannot exclude a 
rise in asymmetry in the Xp < region. The default Pythia/Jetset model, 
which predicts a negative asymmetry over the range > 0, does not provide 
a good description of our data. The data show a trend similar to predictions 
of the two component models MM- 



We have presented data on A^ production asymmetries in both the forward 
{xp > 0) and backward {xp < 0) regions. The range of Xp covered, —0.1 < 
Xp < 0.6 allows the first simultaneous study of the A^ production asymmetry 
in both the negative and positive Xp regions. Our results show a uniform, 
positive asymmetry of (12.7 ± 3.4 ± 1.3) % after acceptance corrections over 
the entire kinematical range studied . 
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Figure 1. pK^TT^ and pi^"'"7r~ invariant-mass distributions for Xp < and .x^ > 0. 
All distributions show a clear A+ and A~ signal. Fitting to a Gaussian signal and 
quadratic background yields 122 zt 17 Aj" (a) and 92 it 15 A~ (b) in the negative Xp 
region and 903 it 42 A+ (c) and 702 it 40 A;: (d) in the positive Xp region. 
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Figure 2. A+ and A;: efficiency ratios, e/e as a function of Xp (upper plot) and 
(lower plot). The efficiency ratios for xp (p|) are integrated over the (xp) range 
of the data set. 
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Figure 3. — A~ asymmetry as a function of Xp (upper plot) and (lower plot). 
Full lines are the prediction of Pythia5.7/Jetset7.4 . The asymmetry for xp (Pt) 
is integrated over the p| (xp) range of the data set. Error bars show the statistical 
and systematic errors added in quadrature. 
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Table 1 

Measured A+ — A~ asymmetry as a fmiction of Xp and pj,. 



Xf 


bin 


Asymmetry 


Pt 


bin 


Asymmetry 


-0.1 


-0.0 


0.20 ±0.10 ±0.06 


0.0 


-1.0 


0.13 ±0.05 ±0.01 


0.0- 


-0.1 


0.11 ±0.05 ±0.03 


1.0- 


-2.0 


0.15 ±0.07± 0.05 


0.1- 


-0.2 


0.14 ±0.05 ±0.03 


2.0 


-3.0 


0.13 ±0.09 ±0.03 


0.2- 


-0.3 


0.14 ±0.10 ±0.05 


3.0 


-5.0 


0.17 ±0.10 ±0.04 


0.3- 


-0.6 


0.30 ±0.39 ±0.15 


5.0- 


-8.0 


0.13 ±0.14 ±0.04 



Table 2 



Comparison of the — A~ asymmetries as measured by the E791, ACCMOR and 
SELEX Collaborations and predicted by three models. Experimental results in the 
Xp > region are for different Xp ranges. See the text for details. 



Xp region E791 


ACCMOR [ 
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SELEX i SF i 


IC i 


2CRi 


x^< 0.20 ±0.10 ±0.06 
Xp> 0.123 ±0.037 ±0.016 


0.005 ±0.079 




+ 

0.25 ±0.15 


+ 



± 
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